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Does the assignment order of a fixed collection of slightly distinct subsystems into given commu¬ 
nication channels influence the overall ensemble behavior? We discuss this question in the context of 
complex networks of non-identical interacting oscillators. Three connection configurations found in 
Nature are considered here: Similar, Dissimilar and Neutral patterns. These strategies correspond 
respectively to oscillators alike, distinct and indifferent relate to its neighbors. To construct such 
scenarios we define a vertex weighted graph measure, the Total Dissonance, which comprises the 
sum of the dissonances between all neighbor oscillators in the network. Our numerical simulations 
show that the more homogeneous is a network, the higher tend to be both the coupling strength 
required to phase-locking and the associated final phase configuration spread over the circle. On 
the other hand, the initial spread of partial synchronization occurs faster for Similar patterns in 
comparison to Dissimilar ones, while neutral patterns are an intermediate situation between both 
extremes. 


Some social and biological studies about multi-agents 
reveal that units tend to select similar peers to inter¬ 
act mm- However, there are systems which behave 
in opposite manner, where their components preferen¬ 
tially choose to connect themselves to others with some 
distinct inner characteristics [3]. In fact, Nature seems 
to favor the former or the later construction, which we 
respectively call Similar or Dissimilar (neighborhood) 
patterns, to achieve different agendas H]. This arti¬ 
cle explores ideas inspired by these scenarios within the 
non-identical phase oscillator Kuramoto model with local 
mean field coupling, which is one of the main paradigms 
to describe collective behavior and synchronization [5]. 
This model is also interesting because it approximates dy¬ 
namics of a large class of non-linear oscillators near limit 
cycle, under weak mutual interaction [6]. Besides, this 
is an active research field with a number of applications 
from different areas CHID], highlighting the fundamental 
role that synchronization plays. 

Our numerical approach is based on a novel vertex 
weighted graph measure: the total dissonance. This 
quantity can be regarded as a generalization of the classi¬ 
cal concept of dissonance, that is, the natural frequency 
difference from the case of two coupled oscillators [B]. 
So, we define Similar and Dissimilar neighborhood pat¬ 
terns as the assignments of non-identical oscillators into 
the coupling graph which yield, respectively, significantly 
lower and higher total dissonance values. Otherwise, if a 
pattern has no strong bias, we call it Neutral. Given a 
fixed choice of inner properties for each oscillator and a 
fixed coupling graph, we search for Similar and Dissimi¬ 
lar patterns via an optimization algorithm interchanging 
oscillator’s positions into the graph nodes. Finally, mas¬ 
sive numerical simulations are performed to grasp the 


influence of these three different neighborhood patterns 
over phase-synchronization quantifiers. Regular, scale- 
free, random, small world and community networks are 
considered m to provide evidences about the ubiquity 
of our argument. 

About related material, Ref. EH and Hi respectively, 
explore first and second order Kuramoto Model versions, 
both including local correlations between oscillator’s nat¬ 
ural frequency and node-degree. They report an explo¬ 
sive synchronization in the first case and cascade syn¬ 
chronization, according to the node-degree, in the sec¬ 
ond. Our methodology introduces a diverse relationship 
between natural frequencies and coupling graph, as will 
be discussed and illustrated in the text. 

Optimization studies also laid foundation to our re¬ 
search. In Ref. H an algorithm is proposed to con¬ 
struct optimized networks related to a combination of 
local and global synchronization measures. Their objec¬ 
tive function is computed and refined after successive nu¬ 
merical integrations. Although we follow a different ap¬ 
proach, we point out that our results also support that 
“the early onset of synchronization and rapid transition 
to the phase-lock are conflicting demands on the network 
topology” M- Ref. |T5] associates a percolation process 
to the spread of synchronization. In addition, they con¬ 
sider node interchange in the graph based on a vertex 
weighted graph measure. However, their characteriza¬ 
tion takes into account only the phase sign of neighbor 
oscillator. Even so, we also found a similar explosive 
synchronization. 

It is a common sense that a way to achieve more ho¬ 
mogeneous neighborhood patterns is to gather members 
with closer intrinsic dynamics into communities. Thus, 
articles that investigate this framework can also benefit 
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from our findings. Ref. 1161 addresses a Kuramoto model 
of identical oscillators showing that, as the transient time 
dies out, synchronization occurs in stages matching the 
granular communities of the coupling graph. Ref. [T7; 
deals with communities of oscillators having essentially 
different natural frequencies. The authors of this paper 
discuss ways to promote or suppress synchrony on indi¬ 
vidual subgroups. Ref. [18] introduces a dynamic feed¬ 
back control to produce intra-communities synchroniza¬ 
tion regarding communities of identical non-linear oscil¬ 
lator. Accordingly, Similar, Neutral and Dissimilar pat¬ 
terns can be a feasible tool to adjust synchronization 
properties within the context of related works. 


MODEL AND METRICS 

We consider a system of N phase oscillators, whose 
dynamics for the *-th oscillator is 


N 



where oj = (uq,..., wjv) € are the oscillator’s natural 
frequencies. We consider ui with zero mean [19j . ran¬ 
domly drawn from the uniform distribution over [— 7r, 7 t] . 
The coupling strength e > 0 is the system parameter that 
adjusts the intensity of attractiveness between neighbor 
oscillators. The symmetrical coupling graph is expressed 
by its adjacency N x N matrix A , so that A,, = 0; 
Aij = 1, if oscillators i, j are neighbors (adjacent); and 
A,j = 0, otherwise. We assume connected graphs, mean¬ 
ing that there is a sequence of edges joining any two 
vertexes in the graph. Also, di := ]Tb v A l:) stands for 
the i-th vertex degree. The Laplacian matrix is defined 
as L := diag(di,..., djv) — A and its eigenvalues are 
0 = Ai<A 2 <...<Ajv. The first non-trivial eigen¬ 
value A 2 , the algebraic connectivity , is greater than zero 
if and only if the graph is connected J2D{. 

On one hand, analytical results ED guarantee con¬ 
vergence to a unique (modulus 27r) stable phase-locked 
regime, where phase differences between every two os¬ 
cillators becomes constant. Precisely, this convergence 
occurs if the coupling strength e is large enough in com¬ 
parison with 11^11 Ajv/A|, where 11.11 denotes the Euclidean 
norm in R w . Because A 2 increases when the graph di¬ 
ameter D is decreased and Ajv decreases with its maxi¬ 
mum degree d max [22], phase-locking can be achieved for 
smaller values of e mostly with the reduction of D , but 
also with smaller values of ||w|| and d max . On the other 
hand, if we consider a system with only two phase oscil¬ 
lators, it is well known [6] that the relation between its 
dissonance v := uq — uq versus the coupling strength e de¬ 
termines the synchronization regimel23j. So, inspired by 
this context, we introduce the total dissonance measure 


for vertex weighted graphs as 
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i,j = 1 


Since we consider symmetrical and connected coupling 
graphs, it is straightforward to check that v To ta.i = 0 
if and only if all oscillator are identical. If we write 
Wotai = J'Totai(w), this measure quantifies how far oj is 
from a condition where all natural frequencies are iden¬ 
tical. Therefore, i/ Tota 1 encompasses information about 
the total spreading of oj by summing up individual dis¬ 
sonances over the coupling graph edges. 

The norm of the global mean field, the order parameter, 
will be denoted by R(9) = 1/N J^iLi el0i ■ This quantity 
R ranges from 0 to 1, respectively indicating that the 
ensemble gradually changes from null global mean field, 
where all phasors e l8i cancel out, to full synchronization, 
where 6\ = ... = On- We also make use of the edge partial 
synchronization index between two oscillators i,j, 
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where t r is a large enough transient time [213 • Oscillators 
i,j are phase-locked, that is 0i(t) — 0j(t) converges to a 
constant value, if and only if Sij = 1. Moreover, if this 
index is decreased towards zero, then weaker forms of 
synchronization and later uncorrelated trajectories occur 
[251 . We average contributions of all neighbor oscillators 
in the network to define the partial synchronization index 
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i,j =1 


(3) 


where E := X3fj=i quantity of directed edges 

in the graph. Of course, the number of undirected edges 
is E := E/2. Thus, 5 = 1 means that the whole ensem¬ 
ble is phase-locked, while S ~ 0 yields very low coherent 
ensemble behavior. Note that R(0(t)) converges to a con¬ 
stant value [26j if and only if S = 1. 

An Adams-Bashforth-Moulton Method for numerical 
integration is applied. Transient time of at least 2.10 3 
units of time was suppressed from the data, while the 
convergence of approximations of S over successive time 
windows of 10 3 units of time was the criterion to inter¬ 
rupt the integration. The mean value of R(0(t)) after 
the transient is denoted by (R). For a given choice of 
parameters and initial conditions, we indicate by e PL the 
smallest critical coupling strength e > 0 inducing phase¬ 
locking, i.e. 5 = 1 and R{0(t)) converges to a constant 
value, which we denote by i? PL - 

Several complex networks topologies ED with N nodes 
and E (undirected) edges are considered: 4-Regular (N 
RE), Barabasi-Albert (N BA), Erdos-Renyi (N ER) and 
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Watts-Strogatz [27] (N WS). Experiments with relatively 
small networks with N = 50 are performed for the sake of 
easy visualization. Larger ones, with N = 500, are also 
addressed to illustrate graphs closer to the theoretical de¬ 
gree distribution El, yet feasible to massive numerical 
integration. To diminish computational cost and to al¬ 
low comparison among network topologies, we consider 
graphs with E = 4 N directed edges, which yields mean 
node-degree (d) = E/N = 4. 

Besides, a real world complex network with community 
structure, denoted by 105 CO, is included in the simu¬ 
lations: the Krebs-Amazon Political Books network 128] , 
This graph comprises N = 105 nodes, E = 882 edges and 
2 communities. 

A single choice of natural frequencies ui £ were ran¬ 
domly drawn, matching network sizes N = 50,105, 500. 
The term Configuration stands for a specific pair of graph 
and permutation of the vector ui. For a fixed coupling 
graph, if one obtains permutations of u> such that the as¬ 
sociated ^Totai is closer to its minimum or maximum, we 
have Similar (S) or Dissimilar (D) patterns, respectively. 
On the other hand, when no optimization process is ap¬ 
plied and z/ Total (w) is far from extremes values, we call it 
a Neutral (N) pattern. 

Fig. 0 illustrates graphs following RE, BA and CO 
topologies with Similar, Neutral and Dissimilar config¬ 
urations. From the RE graph with Similar configura¬ 
tion of this figure, one realizes a gradual transition of uj t . 
It happens because this pattern reduces the dissonance 
through edges, so each node presents natural frequency 
close to the respective average of its neighbors. However, 
this ordering arises differently depending on the network 
topology. Hubs of the BA graph were colored with lighter 
tones, corresponding to the overall mean of natural fre¬ 
quency distribution. But in the CO graph, positive and 
negative natural frequency values were placed into dis¬ 
tinct communities, with hubs close to ±7t/2 and central 
nodes (in between communities) close to null cc,;. 

About the Dissimilar configurations from Fig. [f] the 
opposite organization is found: each node receives natu¬ 
ral frequency far from its neighbors. For the RE network, 
we clearly notice sequences of connected nodes with al¬ 
ternating positive and negatives values of w,;. Besides, 
BA and CO graphs presented connected hubs with larger 
natural frequencies and opposite signs. Eventually, Neu¬ 
tral configuration can be regarded as a blending between 
both previous configurations. 

Let us now describe a mechanism to seek Similar and 
Dissimilar patterns. For a given graph, we use succes¬ 
sive iterations of the Simulated Annealing optimization 
method [25] to search over the AH existing permutations 
of the vector u>. This problem may be challenging by 
itself due to its computational complexity. To enhance 
performance, we betake the observed features described 
in the previous paragraph to empirically tailor our algo¬ 
rithm. 
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FIG. 1. (Color online) Examples of Similar, Neutral and Dis¬ 
similar patterns related to A-nodes graphs following network 
topology: 4-Regular (N RE); Barabasi-Albert ( N BA), both 
with E = 100 edges and; a community graph ( N CO) with 
E = 441. Vertex color according to natural frequency Wi 
ranging from —n (blue) to 7r(red); and vertex size propor¬ 
tional to node-degree. The associate total dissonance i/ Tota i is 
displayed. 


Let a/ 1 ! denote the input to the optimization process. 
This initial seed to search Similar pattern is constructed 
by gathering the values of {u> i,..., lon} which are close 
to each other into (modularity-based) communities 1301 . 
No special node ordering within each community is im¬ 
posed so far, but an auxiliary step of the SA algorithm is 
used to reorder communities to reduce u Total . For Dissim¬ 
ilar patterns, ui 1 ! is allocated by matching larger vertex 
degrees with higher absolute values of{wi,... ,ojn}- 

After that, either to minimize or maximize the total 
dissonance otal , successive iterations of the SA are ap¬ 
plied. For each stage k = 2,..., 30, its initial seed (J k ^ is 
given by the outcome of the previous iteration. Alternat¬ 
ing strategies to randomly evolve approximate solutions 
are employed at each stage. 

First, nodes to be swapped are drawn with probabil¬ 
ity proportional to vertex degree, since they are the ones 
which influence the most ;/ Total . We also pick vertexes to 
be interchanged only within the same modularity-based 
community. Thus, we can focus the optimization process 
on relatively detached graph regions. Lastly, to facilitate 
escape from local critical points and also to fine adjust¬ 
ments, we consider all nodes equiprobable. 
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Patterns: ■ Similar ■ Neutral ■ Dissimilar 


FIG. 2. (Color online) Total Dissonance i/ To tai distribution 
chart of different topologies with Similar (Blue), Neutral 
(Gray) and Dissimilar (Orange) patterns. Data correspond 
to the single members from categories 50 RE, 500 RE and 
105 CO; and 100 different elements for the others. Colored 
bars indicate the overall range of I'Totai obtained, while mean 
values of each distribution are joined by a black line within 
each category. 

Recall that any Configuration (of a graph A and natu¬ 
ral frequencies u>) which suffered no optimization process 
is said to be Neutral. Yet, when we apply to a Configura¬ 
tion (A, ui) the optimization process described above, we 
obtain Similar or Dissimilar patterns with minimization 
or maximization i/ Total , respectively. So, each pair (A, to) 
generates the three patterns. 

Since RE topology is deterministic and the CO graph 
was extracted from a data set, these classes contain a sin¬ 
gle member to be analyzed. Yet, we randomly generate 
and include in our experiments 100 graphs of BA, ER 
and WS network topologies. 

Fig. [2] displays a distribution chart of the total disso¬ 
nances r^Totai obtained for the categories included in this 
article. From this figure, one notices a sharp distinction 
among patterns since there is no iz Total range overlapping 
within each category. 

The set of graphs with N = 500 presented values of 
z'Totai three times smaller than the set with N = 50. Fur¬ 
thermore, both sets were qualitatively alike, which is an 
evidence that we were capable to produce Similar and 
Dissimilar neighborhoods in the large networks, at least 
as well as the small ones. 

Because there is only a normalization by the number 
of vertex in the graph, but not by the quantity of edges, 
in Eq. ([2]), the total dissonance decreases with smaller 
mean degrees in all categories as expected. [31] BA, ER 
and WS topologies were almost indistinguishable in the 
Neutral pattern, but were slightly higher in this order, for 
the Dissimilar case. All patterns of 50 RE and 105 CO 
graphs yielded smaller but comparable with each other 
total dissonance values. 

We remark that the essential purpose of the total dis¬ 
sonance Mpotai i n this work was to lay foundation for the 
definition and search of Similar, Neutral and Dissimilar 


configurations. We could have formulated the subjective 
concepts of Similar (homogeneous) and Dissimilar (het¬ 
erogeneous) networks based on variations of Eq. ([2]) . It 
is deferred for awhile to future research the study of cor¬ 
relations between gradual transitions of the values of this 
metric, far from the Similar and Dissimilar extremes, and 
synchronization features. 

We abuse notation and also denote by S, (R), e PL and 
i? PL the associated mean values of these synchronization 
quantifiers considering all graphs of each category. A 
fixed random choice of initial condition 9° £ R N is drawn 
from a uniform distribution over the unit circle for each 
network size N. So, Fig. [3] displays the values of S 
and ( R ), as solid and dashed line resp., obtained through 
numerical integration for the three neighborhood config¬ 
urations colored like in Fig. [2] The time variable of both 
S and (R) times series are rescaled to end at the associ¬ 
ated mean critical coupling value er PL . The S lines finish 
at value 1, within the numerical tolerance, while the fi¬ 
nal value of (R) lines equal to the mean critical order 
parameter i? PL . 

First, one focus on the phase-locking measures e PL and 
i? PL . Irrespective to network size, overall results were 
alike. In general, e PL decreases from Similar, Neutral to 
Dissimilar patterns; while f? PL tend to increase in the 
same ordering. In all topologies, Similar cases demanded 
higher coupling strength to achieve phase-locking. In 
particular since this holds true even for RE networks, 
it shows that total dissonance patterns induce a different 
phenomena than the ones from Refs [T2J and [T51 

Moreover, even when these networks phase-locked, R 
converged to smaller values of f? PL . In other words, Sim¬ 
ilar ensembles tend to be harder to synchronize and to 
converge to regimes where oscillators were more spread 
around the unit circle than their counterparts. Neutral 
patterns required smaller £ PL than Dissimilar ones. RE 
graphs were the only exception for these behavior of f? PL . 

For all topologies, higher values of e PL were measured 
when N was multiplied by 10. On the other hand, 
larger networks yielded higher f? PL for Similar and Neu¬ 
tral neighborhoods, but slightly smaller f? PL for Dissimi¬ 
lar ones. 

At this point, the influence of Similar, Neutral 
and Dissimilar patterns over the emergence of phase- 
synchronization is investigated, specially related to cou¬ 
pling strengths e much smaller than e PL . 

Again, except for RE graphs, we verify that Similar 
patterns favor weaker synchronization regimes, since the 
initial growth of S and (R) for small coupling strength e 
is more prominent. However, beyond intermediate val¬ 
ues of c, Dissimilar patterns surpass the Similar ones 
through an abrupt transition. The Neutral case is in- 
between these two extremes, closer to the behavior of 
the Dissimilar group. If we compare network topologies, 
BA and ER graphs displayed close values of S, which 
were smaller than WS ones for small and intermediate 
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<R>: Similar Neutral Dissimilar 

S: - Similar Neutral Dissimilar 


FIG. 3. (Color online) Mean order parameter after transient (R) and partial synchronization index S, solid and dashed lines 
resp., as a function of the coupling strength e for different graph topologies. Average values of all graphs simulated within each 
category are shown. Similar, Neutral and Dissimilar cases are respectively plotted in Blue, Gray and Orange. Lines are drawn 
up to e equal to the respective average critical phase-locking e PL . 


values of e. 

A parallel of our findings could be made with percola¬ 
tion of conflicting ideas, associating communication and 
agreement with the emergence of synchronization and 
phase-locking, respectively. In Similar scenarios, inter¬ 
action mostly occurs among people with closely related 
culture backgrounds. Thus, communication can easily 
spread locally, but the overall population, which contains 
diverse members, will hardly find a compromise. On the 
other hand, when networks contain more heterogeneous 
neighbors, as in the Neutral and Dissimilar cases, com¬ 
munication demands higher effort to be established. But 
after that, the whole ensemble is capable to rapidly reach 
consensus. 

In summary, experiments with several network topolo¬ 
gies were analyzed and strong numerical trend was found. 
The Neutral case behaves in general between both ex¬ 
tremes, closer to the Dissimilar case. Except for RE 
networks, under small coupling strength e, Similar pat¬ 
terns yield larger values of partial synchronization index 
S, meaning early synchronization ongoing. In contradis¬ 
tinction, Dissimilar ones present smaller values of S, but 
undergo abrupt increment until phase-locking. Moreover, 
all networks with Similar patterns required higher val¬ 
ues of coupling strength to achieve phase-locking, while 
Dissimilar patterns converged to regimes closer to full 
synchronization. 
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